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We introduce a parametrization scheme for J/tp(tp') — > VP where the effects of SU(3) flavor sym- 
metry breaking and doubly OZl-rule violation (DOZ1) can be parametrized by certain parameters 
with explicit physical interpretations. This scheme can be used to clarify the glueball-gq mixing 
within the pseudoscalar mesons. We also include the contributions from the electromagnetic (EM) 
decays of J/tp and tp' via J/ip(ip') — > 7* — > VP. Via study of the isospin violated channels, such 
as J/ip(ip') — + pn, prj ', lott and 4 >7T °i reasonable constraints on the EM decay contributions are 
obtained. With the up-to-date experimental data for J/ip(tp') — + VP, J/xp(ip') — » 7P and P — > 77, 
etc, we arrive at a consistent description of the mentioned processes with a minimal set of parame- 
ters. As a consequence, we find that there exists an overall suppression of the ip' — > 3g form factors, 
which sheds some light on the long-standing u pn puzzle" . By determining the glueball components 
inside the pseudoscalar n and n' in three different glueball-gg mixing schemes, we deduce that the 
lowest pseudoscalar glueball, if exists, has rather small qq component, and it makes the ??(1405) a 
preferable candidate for glueball. 



PACS numbers: 13.20.Gd, 12.40.Vv, 13.25.-k, 12.39.Mk 



I. INTRODUCTION 



Charmonium decays into light hadrons provides unique places to probe light hadron structures. In 
particular, in the hadronic decays of charmonia such as J /ip, rj c and Xcj etc, the annihilation of the heavy 
cc pair into intermediate gluons, which must be then hadronized into hadrons, could favor the production 
of unconventional hadrons such as glueball and hybrids in the final-state. Such states, different from the 
conventional qq or qqq structures in the non-relativistic constituent quark model, can serve as a direct 
test of QCD as a non-Abelian gauge theory. 

In the past decades, the exclusive hadronic decays of J/ip(ip') — > VP have attracted a lot of attention 
in both experiment and theory. These processes, in which the helicity is non-conserved, must be sup- 
pressed according to the selection rule of pQCD hadronic helicity conservation due to the vector nature 
of gluons [l[. According to the "pQCD power" suppression of helicity conservation and breaking, one 
should have BR(tp' — ► pir) / 'B R(J / 'ip — > pn) ~ (Mj/^/M^>) 6 ~ 0.35 l]. However, the experimental data 
show that this ratio is badly violated in reality, BR{ip' — > pn)/BR(J/ip — * pn) ~ (0.2 ± 0.1)% 0, i.e. 
much greatly suppressed. It led to the so-called "pn puzzle" in the study of J I ip and j>' exclusive decays, 
and initiated a lot of interests to the relevant issues 0, H i, H, S II E E3, Ql E 111 0, E GJ, G3- 
An alternative expression of the " pn puzzle" is via the ratios between J jvb and if)' annihilating into three 
gluons and a single photon: 

_ BR(ijj' — > hadrons) 
R = 



BR(J/ip — > hadrons) 
BR(ip' -> e+e~) 



BR { J,^^ l2% ' W 
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which is empirically called "12% rule". Since most of those exclusive decays for J/ip and ip' seem to 
abide by this empirical rule reasonably well, it is puzzling that the ratios for pir and K*K + c.c. deviate 
dramatically from it. A recent article by Mo, Yuan and Wang provides a detailed review of the present 
available explanations (see Ref. [l8| and references therein). 

A catch- up of this subject is an analysis by authors here for the electromagnetic decays of J/ip(ip') — > 
7* — > VP [IjSj]- There, it is shown that although the EM contributions to J/ip — > VP are generally small 
relative to the strong decays as found by many other studies [201 ] , they may turn to be more competitive 
in tp' — > VP compared with ip 1 — > 3g — > VP amplitudes, and produce crucial interferences For a 

better understanding of the " pn puzzle", a thorough study of J/tp(ip') — > VP including both strong and 
EM transitions and accommodating the up-to-date experimental information [2l|, H2, l23l . l2~j ] should be 
necessary. 

The reaction channels J/t/j(ijS) — » VP also give access to another interesting issue in non-perturbative 
QCD. They may be used to probe the structure of the isoscalar '", i.e. rj and 7/ in their recoiling isoscalar 
vector mesons u> or (f>. Since vector meson u and (j> are almost ideally mixed, i.e. u> = (uu + dd)/^/2 and 
4> = ss, the decay of J/i/KV'') ~~ * wrj' \ ( fi r l> and 4>r]' will provide information about 77 and rj' , which can 
be produced via the so-called singly OZI disconnected (SOZI) processes for the same flavor components, 
and doubly OZI disconnected (DOZI) processes for different flavor components. 

Nonetheless, it can also put experimental constraints on the octet-singlet mixing as a consequence of 
the U(1) A anomaly ofQCD [H, H3, [H, [H, [H OH [H, [H . Recently, Leutwyler (j^, and Feldmann 
and his collaborators [331 ] propose a 77-77' mixing scheme in the quark flavor basis with an assumption that 
the decay constants follow the pattern of particle state mixing, and thus are controlled by specific Fock 
state wavefunctions at zero separation of the quarks, while the state mixing is referred to the mixing in 
the overall wavefunctions. In the quark flavor basis it is shown that higher Fock states (specifically, \gg)) 
due to the anomaly give rise to a relation between the mixing angle and decay constants. This analysis 
brings the ques tion of gluon components inside rj and 77' to much closer attention from both experiment 
and theory [33l l34j. and initiates a lot of interests in the pseudoscalar sector, in particular, in line with 
the search for pseudoscalar glueball candidates. 

The lattice QCD (LQCD) calculations predict a mass for the lowest pseudoscalar glueball around 2.5 
GeV 35j , which is higher than a number of h resonances observed in the range of 1~ 2.3 GeV. However, 
since the present LQCD calculations are based on quenched approximation, the glueba ll sp ectrum is still 
an open question in theory. In contrast, QCD phenomenological studies [36|, [33, [H, [39| favor a much 
lower h glueball mass, and make r/(1405) a good glueball candidate due to its strong couplings to KKtt 
and ao(980)7r and absence in 77 — > KKn and rjinr [2J. Study of J/tp — > VP in flavor parametrization 
schemes was pursued in Ref. [20j , where the rj and rj 1 were treated as eigenstates of qq. By assuming the 
77(1405) to be a glueball essentially and introducing a DOZI suppression factor for </>G and ujG amplitudes 
relative to the SOZI processes, the authors seemed to have underestimated the branching ratios for 
BR(J/tp -> 077(1405)) and BR(J/ip -> wry(1405)). In Ref. [H a mixing scheme for 77, 77' and 7/(1405) 
in a basis of (uu + dd)/^/2, ss and G (glueball) was proposed and studied in J/ip — > VP. 7/(1405) was 
found to be dominated by the glueball component, while the glueball component in 77' was also sizeable. 
In both studies, the EM contributions were included as a free parameter in J/ip —* VP. However, as we 
now know that the EM contributions in are important [Hj]. Thus, a coherent study of J/ip(^') — > VP 
including the EM contributions would be ideal for probing the structure of the pseudoscalar mesons. 

On the theoretical side there is only one possible term in the effective Lagrangian for J/ip(ip') — > VP 
and it will reduce the number of the parameters needed in phenomenological studies. On the experimental 
side many new data with high accuracy are available. Hence, in this work, we shall revisit J /^(tp 1 ) — * VP 
trying to clarify the following points: (i) What is the role played by the EM decay transitions? By isolating 
and constraining the EM transitions in a VMD model, we propose a parametrization scheme for the strong 
transitions where the SU(3) flavor symmetry breaking and DOZI violation effects can be accommodated. 
A reliable calculation of the EM transitions in turn will put a reasonable constraint on the strong decay 
transitions in J/ip(ip') — > VP, and thus mechanisms leading to the "pir puzzle" can be highlighted, (ii) 
We shall probe the glueball components within 77 and 77' based on the available experimental data and 
different quarkonia-glueball mixing schemes, and predict the glueball production rate in J/ip(ip') 
hadronic decays. 

As follows, we first present the parametrization scheme for J/ip(tp') — > VP and introduce the glueball 
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components into the pseudoscalars. We then briefly discuss the VMD model for J/ip(ip') — > 7* — > VP. 
In Section III, we present the model calculation results with detailed discussions. A summary will be 
given in Section IV. 

II. MODEL FOR J/il>(i>') -» VP 

In this Section, we first introduce a simple rule to parametrize out the strong decays of J/i(>(ip') — > VP. 
We then introduce a VMD model for the J/ip(ifj') EM decays into VP. A recent study of J/ip(tp') — ► 
7* — ► VP shows that the EM decay contributions become important in ip' decays due to its large partial 
decay width to e + e~ though their importance is not so significant in J/ip — > VP [lj|. In the isoscalar 
pseudoscalar sector, we introduce three different glueball-qq mixing schemes for 77, rf and a glueball 
candidate 77". Its production in J/tfj(ip r ) hadronic decays can then be factorized out and estimated. 

A. Decay of J/tp(ip') — ► VP via strong interaction 

The strong decay of J/ip(ip') —f VP via 3g can be factorized out in a way similar to Ref. (4lj . 
First, we define the strength of the non-strange singly OZI disconnected (SOZI) process as 

9J/W) = {(qq)v(qq)p\V \J/i>(^)), (2) 

where Vo denotes the 3<? decay potential of the charmonia into two non-strange qq pairs of vector and 
pseudoscalar via SOZI processes. But it should be noted that the subscription V and P do not mean 
that the quark-antiquark pairs are flavor eigenstates of vector and pseudoscalar mesons. The amplitude 
9j/4>(4>') is proportional to the charmonium wavefunctions at origin. Thus, it may have different values 
for J /if) and ip'. 

In order to include the SU(3) flavor symmetry breaking effects in the transition, we introduce 

R = ((qs) v (sq)p\V \J/ip(i}'))/gj/^>) = ((sq)v(qs)p\V)\ J/ip(i>'))/gj/^>), (3) 

which implies the occurrence of the SU(3) flavour symmetry breaking at each vertex where a pair of ss is 
produced, and R = 1 is in the SU(3) flavour symmetry limit. For the production of two ss pairs via the 
SOZI potential, the recognition of the SU(3) flavor symmetry breaking in the transition is accordingly 

R 2 ~ ((ss) v (ss)p\Vo\J/^'))/g J/ ^ n . (4) 

Similar to Ref. [4l[, the DOZI process is distinguished from the SOZI ones by the gluon counting rule. 
A parameter r is introduced to describe the relative strength between the DOZI and SOZI transition 
amplitudes: 

r = ((ss) v (qq)p\V 1 \J/tp(ip'))/g J / 4 ,^ l) = {(qq)v(ss) P \Vi\ J/ip(il>'))/gj/W), (5) 

where V\ denotes the charmonium decay potential via the DOZI processes; In the circumstance that the 
OZI rule is respected, one expects \r\ ~ 0. 

Through the above definitions, we express the amplitudes for J/tpii/j') — * V(qq)p as 



M K , +K - 


- (K*+K-\V \J/il>W)) =9J/W)R , 


(6) 


M p + W - 


= (p + ir~\V \J/ip(il)')) =5j/v(V) ' 


(7) 


M p0jr o 


= (p°AVo\J/W)) =9j/w) , 


(8) 




= (4>{ss) P \Vo + V 1 \J/^'))=R 2 g.j / ^, ) {l + r), 


(9) 


M^nn) 


- Wnn)p|Vi|J/V»(^)) = y/2r{4>{aa)p\Vo\JH>W)) = y/2rRg mm , 


(10) 




- {u(8s) P \Vi\J/il>{%l>')) - y/Wrgj/w), 


(11) 




= (u(nn) P \Vo + Vi|JMV0) = 9J/^)0- + 2r )> 


(12) 
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where the channels KK* and p 7r + have the same expressions as their conjugate channels; nn = (uu + 
dd)/V2 is the non-strange isospin singlet. We have assumed that <\> and w are ideally mixed and they are 
pure ss and nn, respectively. 

The recoiled pseudoscalars by <fr and lu can be either n or n' for which we consider that a small glueball 
component will mix with the dominant qq in the quark flavour basis. The detailed discussion about the 
qq and glueball mixings in n and rf will be given in the next section. Here, we relate the production of 
glueball G in J /ip{ip') — > loG and 4>G with the basic amplitude gj/^Uj,') by assuming the validity of gluon 
counting rule in the transitions [ill. |42|, i.e. 

((qq) v G\V 3 \ J/V#')> = 9j/W) , (13) 

where V3 denotes the glueball G production potential recoiling a flavour singlet qq. This can be regarded 
as reasonable since generally the glueball does not pay a price for its couplings to gluons. Consequently, 
we have 

AUg = (d>G\V 3 \J/^')) = g Jfi>m R 

M uG = {wG\V 3 \J/rl>W))=V2gj/w)> ( 14 ) 

with the SU(3) flavor symmetry breaking considered for the 4> production. 

The above parametrization highlights several interesting features in those decay channels. It shows 
that the decay of J/ip(tp') — ► pn is free of interferences from the DOZI processes and possible SU(3) 
flavor symmetry breaking. Ideally, such a process will be useful for us to determine gj/^(^>). The decay 
of J/ifj(ip') — > K*K + c.c. is also free of DOZI interferences, but correlates with the SU(3) breaking. These 
two sets of decay channels will, in principle, allow us to determine the basic decay amplitude gj/^u,<) 
and the SU(3) flavor breaking effects, which can then be tested in other channels. For 4>n, <pn' , ton and 
tun' , the transition amplitudes will also depend on the 77-77' mixing angles and we will present detailed 
discussions in later part. 

In the calculation of the partial decay width, we must apply the commonly used form factor 

T 2 (P) = |P| 2i exp(-P 2 /8/3 2 ), (15) 

where P and the I are the three momentum and the relative orbit angular momentum of the final- 
state mesons, respectively, in the J/ip(ip') rest frame. We adopt (3 — 0.5GeV, which is the same as 
Refs. [42], [H, H ]. At leading order the decays of J/ip(ip') — > VP are via P-wave, i.e. 1 = 1. This form 
factor accounts for the size effects from the spatial wavefunctions of the initial and final-state mesons. 



B. Decay of J/tp(^') -» VP via EM interaction 

Detailed study of J/^(ip') —> 7* — > VP in a VMD model is presented in Ref. [19(. In this process, 
three independent transitions contribute to the total amplitude as shown by Fig. [21 The advantage of 
treating this process in the VMD model is to benefit from the available experimental information for all 
those coupling vertices. As a result the only parameter present in the EM transition amplitude is the 
form factor for the virtual photon couplings to the initial (J/ip and ip') and final state vector mesons (ui, 
(f>, p and K*). As shown in Ref. [19(, the isospin violated channels, pn, pn' , lutt and 07r, provides a good 
constraint on the 7* VP form factors without interferences from the strong decays. Therefore, a reliable 
estimate of the EM transitions can be reached. 

Following the analysis of Ref. [H[, the invariant transition amplitude for J/ip — * 7* — > VP at tree level 
can be expressed as: 

Mem = Ma+Mb +M C 

€ 9V1 ^a + -^^n + T-^TF 2 ^) e^PVrirvep, (16) 



fv'2 My I fvi Mv2 fvifv2 Mp 
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where e//y denotes the 7*^ couplings, and <?y 7 p is the coupling determined in the radiative decay of 
V — ► 7V; !F a , jFf, and .F c denote the form factors corresponding to the transitions of Fig. For Fig.[2^a) 
and (b) a typical monopole (MP) form factor is adopted: 

= (17) 

with A = 0.542 ± 0.008 GeV and A = 0.577 ± 0.011 GeV determined with a constructive (MP-C) or 
destructive phase (MP-D) between (a) and (b), respectively. We think that the non-perturbative QCD 
effects might have played a role in the transitions at J/ip energy. For instance, in Fig. Eta) and (b) a 
pair of quarks may be created from vacuum as described by 3 Pq model, the pQCD hadronic helicity- 
conservation due to the vector nature of gluon is violated quite strongly. A monopole-like (MP) form 
factor should be appropriate for coping with the suppression effects, and it is consistent with the VMD 
framework. In principle such a form factor can be tested experimentally via measuring the couplings of 
the processes J/ip(if)') — ► Pe + e~ and e + e _ — > -Pe + e~~, respectively, when the integrated luminosity at 
J/ip and the suitable energies for e + e~ colliders is accumulated enough. 

Due to an error of missing a factor of y/2 in Tab. I of Ref. [l9[, we give here in Table Q] the coupling 
e/fv again. Also, we clarify that although the overall factor of \f2 will change the fitted cut-off energy 
A and the fitted branching ratios slightly in Ref. [19], the pattern obtained there retains and the major 
conclusion is intact. 

The form factor for (c) has a form of: 

Tc{ £^ = (l-g2/A2)(l- 9 2/A2) ' (18) 

where q\ = M V1 and g| = My 2 are the squared four-momenta carried by the time-like photons. We 
assume that the A is the same as in Eq. (|17p . thus, T c = T a J~b- 

A parameter 6 is introduced to take into account the relative phase between the EM and strong 
transitions: 

M=M 3g + e lS M EM ■ (19) 

It will then be determined by the experimental data in the numerical fitting. In the limit of 8 = 0, the 
relative phase reduces to the same ones given in Ref. [2fjj . 



C. Mixing of the 77 — 77' and implication of a + glueball 

The 77— rj mixing is a long-standing question in the literature. Here we would like to study the "mixing" 
problem by including empirically a possible glueball component in the 77 and 7/ wavefunctions. We extend 
the mixing of the the 77 and 77' as a consequence of the flavor singlet nn, ss and glueball mixing. The 
corresponding glueball candidate is denoted as 77". In the quark flavor basis, treating 77, 77' and 77" as the 
eigenstates of the mass matrix M with the eigenvalues of their masses M v , M v r and M v 'i, we have 



1 f M * 








UMU' 1 = 


M v , 





V 








(20) 



where U is the state mixing matrix: 





/ 77,77, 




= U I ss 


r\") 


\ G 




(21) 



Three mixing schemes are applied to determine the mixing matrix elememts. 
I) CKM approach 



G 



By assuming that the nh, ss and glueball G make a complete set of eigenstates, and physical states 
are their linear combinations, we can express the mixing in the same way as the CKM matrix with the 
phase (5 = (no CP violation is involved): 

/ C12C13 S12C13 S13 \ 

U = I -S13C23 - C12S23S13 C12C23 - S12S23S13 S23C13 J , (22) 
V S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 / 

where Cij = cos 8ij and Sy = sin 9ij with 8ij the mixing angles to be determined by the experimental data. 
The feature of this approach is that the completeness guarantees the unitary and orthogonal relations 
of the matrix. However, it also implies that mixings beyond qq and glueball are not allowed. This may 
be a strong assumption since resonances and exotic states such as tetraquarks with the same quantum 
number could also mix with the qq and glueball. Because of this, the CKM approach will test the extreme 
condition that only ground state qq and G mix with each other. 
II) qq-G mixing due to higher Fock state \gg) 

This scenario is initiated by the QCD U(1)a anomaly. In a series of studies by Feldmann et al. [H, HH, 
it is pointed out that an appropriate treatment of the mixing requires to distinguish matrix elements 
of 77 and 77' with local currents and overall state mixings. Nonetheless, the divergences of axial-vector 
currents including the axial anomaly connect the short-distance properties, i.e. decay constants, with 
the long-distance phenomena such as the mass- mixing, and highlight the twist-4 \gg) component in the 
Fock state decomposition. In the quark flavor basis, higher Fock state due to anomaly can give rise to a 
non- vanishing qq — > \gg) transition. Similar to the prescription of Ref. [HI, l34j. we introduce the glueball 
components in r\ and n' via the higher Fock state decompositions in rj n and n s : 

\r) n ) = V q \nn) + V g g \gg) + ---, 

\ Vs ) = * s \ss) + W a \gg) + ---, (23) 

where ^ q and \P S are amplitudes of the corresponding qq while ty a q and are those of the \gg) components; 
The dots denote higher Fock states with additional gluon and/or qq which is beyond the applicability of 
this approach. The presence of the gluonic Fock states in association with isoscalar nn and ss breaks the 
orthogonality between nn and ss, and we parametrize such an effect in n and n' wavefunctions: 

I?/) = [a(cos <fi\nn) — sin</>|ss)) + 6(cos0 — smcj))\gg)}, 



\n') = {a(sm(f)\nn} + cos0|ss)) + 6(sin0 + cos (j))\gg)], (24) 

'No 



where the normalization factors are Ni = a 2 + b 2 (l — sm2(p) and N2 = a 2 + 6 2 (l + sin20), and parameters 
a and b are to be determined by experiment. In this treatment, a and b now correlate with the mixing 
angle = + arctan \^2, with 9 as the octet-singlet mixing angle defined in the SU(3) symmetry limit 
(6 — > 0). For the commonly accepted range 9 ~ —24.6° or —11.5° from the linear or quadratic mass 
formulae, respectively the glueball component in n has a strength of 6(cos(/> — sin <fi)/ y^Vi , which is 
relatively suppressed in comparison with that in 7/, i.e. b(smcj) + cos <fi) / V ^2 , for 0° < (j> < 90°. This 
naturally gives rise to the scenario addressed in Ref. [H, Hj] . 

Based on the unitary and orthogonal relation for three mixed states, we can derive: 

= i-{ x j + 4), 

= i-{y\ + yl), (25) 
= 1 -(£ + 4)- 

which will provide information about the qq and glueball components in n' though it is not necessary for 
the unitary and orthogonal relation being satisfied. If other higher Fock states which are negligible in r\ 
and 77' are present in n" with sizable amplitude, further constraints on the mixing wavefunction for 77" 
will be needed. 

Ill) Mixing in an old perturbation theory 
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Considering the mixing between quarkonia and glueball at leading order of a perturbation potential, 
i.e. the transition strength between two states are much smaller than the mass difference of these two 
states, we can then express the physical states as [HI 



\V) = 



1 

C~i 



V2f b 



V2f a 



l 

C T 2 



M nn - M sg 
V2fb 



M nn - M G 



\G) 



W) = TT 



I 

c 3 



M,g - Mr, 



M n - Mr, 



nn) + \ss)+ \G) 



M sS - Mq 



Mq - M,. 



\G) 



(26) 



where f a = (ss\V g \G) — (nn\V g \G) / V 2 is the mixing strength for glueball-gg transitions, while 
ft = (qq\V q \ss) is the ss and non-strange qq mixing strength via transition potential V q . 61,2,3 are 
the normalization factors. M n n, M sS and Mq are masses for the pure h quarkonia and glueball states, 
respectively. They will be determined with parameters /„ and fb by fitting the experimental data to 
satisfy the the physical masses for 77, 77' and 77", and unitary and orthogonal relations. 

The above three parametrization schemes address different aspects correlated with the glueball-gg 
mixing. The CKM approach automatically satisfies the unitary and orthogonal relations, and allows the 
mixing exclusively among those three states. In reality, this may not be the case, and other configurations 
with the same quantum number may also be present in the wavefunctions. The second scheme introducing 
quarkonia-glueball mixing within 77 and 77' via higher Fock state. In principle, there is no constraint on the 
77" configuration. Thus, unitary and orthogonal relations do not necessarily apply to those three states. 
The third scheme applies the old perturbation theory and considers the non- vanishing transitions between 
those pure states. As a result, a physical state will be a mixture of quarkonia-glueball configurations. The 
unitary and orthogonal relations will be a constraint in the determination of the parameters. Comparing 
these three different parametrization schemes, we expect that the quarkonia-glueball mixing mechanism 
and implication of the pseudoscalar glueball candidate can be highlighted. 

In Table HT1 the transition amplitudes for the strong decays of J/ip(ip r ) —> VP are given. 



III. NUMERICAL RESULTS 



With the above preparations, we do the numerical calculations and present the results in this section. 



A. Parameters and fitting scheme 



The parameters introduced in these three different approaches can be classified into two classes. One 
consists of parameters which are commonly defined in all three schemes, such as gj/^upi), r and R, 
and S. Parameter gj/^upi) is the basic transition strength for J/ip(ip') — > VP, and proportional to the 
wavefunctions at origin for J/ipty'). Obviously, it has different values for J ftp and ip' decays, respectively. 
Parameter R and r are the relative strengths of the SU(3) flavor symmetry breaking and DOZI violation 
processes to gj/^N,')- They can also have different values in J/ip and ip' decays. Parameter 8 indicates 
the relative phases between the strong and EM transition amplitudes. 

The other class consists of parameters for quarkonia-glueball mixings, for instance, the mixing angles 
in the CKM approach, and masses M nfL , M sS , and Mq for the v quarkonia and glueball, respectively. 
These parameters depend on the mixing schemes, and as mentioned earlier they give rise to different 
scenarios concerning the 77 and 77' mixings. We shall discuss their properties in association with the 
numerical results in the next section. 

It should be noted that in the VMD model for the EM decay transitions [lj| all couplings are deter- 
mined independently by accommodating the experimental information for V\{Vi) — > e + e~ , P — > 77, and 
Vi(V^) — > 7-P or?-^ 7V2, where V\ = J/^p or ifj' , V2 = <i>, p, or K* , and P = n, 77, 7/ or K. This 
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is essential for accounting for the EM contributions properly. Meanwhile, the radiative decays, such as 
V —> 7?7 and 777', rj — > 77 and rj 77, can probe the qq structure of the vector and pseudoscalar mesons. 
Their constraints on the rj and rj mixing have been embedded in constraining the EM transitions. As 
shown in Ref. [l9[, all the experimental data for Vi(I4) — * e + e~, P — > 77, and Vi(V2) — * jP or P — > 7V2 
have been included. 

Given a reliable description of the EM transitions, we can then proceed to determine the strong transi- 
tion parameters and configuration mixings within the pseudoscalars by fitting the data for J/ip(ip') — > VP. 
The numerical study step is taken as follows: i) Fit the data for J/ip(ip') — > VP with only the strong 
decay transitions; ii) Fit the data for J/ip(ip') — > VP including the EM decay contributions. Comparing 
those two situations, information about the role played by the EM transitions, and their correlations 
with the strong decay amplitudes can thus be extracted. The rj and rj configurations can also be con- 
strained. We emphasize that a reasonable estimate of the EM contributions is a prerequisite for a better 
understanding of the underlying mechanisms in J/ip(ip') — > VP. 

B. Numerical results and analysis 

For each of these three parametrization schemes, three cases are examined: i) with exclusive contri- 
butions from J/ip(ip') strong decays; ii) with a MP-D form factor for the EM transitions; and iii) with 
a MP-C form factor for the EM transitions. The parameters are listed in Table IHUfVl and the fitting 
results are listed in Tables IVlTl and [Villi for J/ip and ip' , respectively. 

In general, without the EM contributions the fitting results have a relatively larger \ 2 value. With the 
EM contributions, the results are much improved for both MP-D and MP-C model. To be more specific, 
we first make an analysis of the commonly defined parameters, i.e. gj/^(^>), r, R and 5, and then discuss 
the quarkonia-glueball mixing parameters for each scheme. 

For those commonly defined parameters their fitted values turn to be consistent with each other in those 
three schemes. It is interesting to compare the fitted values for the basic transition amplitude gj/^^') 
for J/ip and ip' . It shows that the inclusion of the EM contributions will bring significant changes to this 
quantity in ip' decays, while it keeps rather stable in J/ip decays. This agrees with the observation that 
the EM contributions in J/ip hadronic decays are less significant relative to the strong ones [20, 40]. The 
absolute value of gj/^(^>) for tp' is naturally smaller than that for J/ip. Since gj/ip(ip') is proportional to 
the wavefunctions at origins [19(, the small fraction, [g^> /gj/^] 2 — 0.018, implies an overall suppression 
of the ip' — > 3g form factors in ip' — > VP, which is smaller than the pQCD expectation, ~ 12%. It is 
worth noting that this suppression occurs not only to pir, but also to all the other VP channels. We shall 
come back to this point in the later part. 

Parameter r, denoting the OZI-rule violation effects, turns to be sizeable in J /if) decays but smaller in 
ip' decays though relatively large uncertainties are accompanying. This is consistent with our expectation 
that the DOZI processes in ip' decays will be relatively suppressed in comparison with those in J/tp 
decays. However, it is noticed that r has quite large uncertainties in tp' decays though the central values 
are small. This reflects that data for ip' decays still possess relatively large errors, and increased statistics 
may better constrain this parameter. 

The SU(3) breaking parameter exhibits an overall consistency in the fittings. Relatively large SU(3) 
flavor symmetry breaking turns to occur in the J/ip decays compared with that in ip' . In the case without 
EM contributions, the SU(3) breaking in ip' decays also turns to be large. This may reflect the necessity 
of including the EM contributions. It shows that the SU(3) symmetry breaking in J/ip can be as large 
as about 34%, while it is about 5-20% in ip' decays. 

The phase angle S is fitted for J/ip and ip' , respectively, when the EM contributions are included. It 
shows that J/tp favors a complex amplitude introduced by the EM transitions. In contrast, the strong 
and EM amplitudes in ip' decays is approximately out of phase. This means that the strong and EM 
transitions will have destructive cancellations in pir, K* + K~ +c.c, but constructively interfere with each 
other in K*°K° + c.c. As qualitatively discussed in Ref. [l9j], this phase can_lead to further suppression 
to ip' — > pir, and also explain the difference between K* + K~ + c.c. and K*°K° + c.c. The relative phases 
are in agreement with the results of Ref. [13] , but different from those in Ref. [IB] • 
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1. CKM approach 

In this mixing scheme it shows that mixing angles 612 and #23 are well constrained by the experimental 
data while large uncertainties occur to #13 when the EM transitions are included. The corresponding 
matrix element sin 612 favors a small value which implies small glueball component inside 77 meson. The 
mixed wavefunctions obtained with different EM interferences are: 

( \rj) = 0.915|rm) ™ 0.404|ss) + 0.003|G) 
without EM : I |tj') = -0.398|nn) - 0.903|ss) - 0.163|G) ; (27) 
[\t]") = O.Q68\nn) + 0.148|ss) - 0.987|G) 



1 77) = 0.904|nrl) - 0.427|ss) + 1.44 x 10~ 5 |G) 
MP-D Model: { \if) = -0.421 |nri) - 0.892|ss) - 0.166(G) ; (28) 

IV) = 0.071 1 wi) + 0.150|ss) -0.986|G) 



and 



\r)) = 0.901|nn) - 0.433|ss) + 5.0 x 10- 7 |G) 
MP-C Model: { \rf) = -0A27\nn) - 0.888|ss) - 0.168|G) . (29) 
\rp) = 0.073|nn) +0.151|ss) - 0.986|G) 

It shows that though 813 has large uncertainties, its central values indicate a small glueball component 
in ?/. In particular, we find sin $13 ~ in the MP-D model, while sin 6*13 = 0.003 or 5.0 x 10~ 7 from the 
calculations without EM contributions or in the MP-C model, respectively. In contrast, the amplitude of 
the glueball component in the rj' wavefunction is much larger. In all three models (and also in all three 
mixing schemes), a stable glueball mixing magnitude of ~ 17% is favored. 

The automatically satisfied unitary and orthogonal conditions lead to the prediction of the structure 
of the i]" . Assuming 77 (1405) corresponding to r)", the CKM scheme leads to a prediction of glueball- 
dominance inside 77(1405). About 15% of ss and 7% of nn are also required in 77" and they favor to be 
out of phase to the G component. 



2. Quarkonia- glueball mixing due to higher Fock state 

Table IIVI shows that the octet-singlet mixing angle 9 for 77 and rf is within the reasonable range 
of —24.6° ~ —11.5° as found by other studies 0]. Both 77 and 77' can accommodate a small glueball 
component in association with the dominant qq. Parameter a and b are fitted by quite different values 
with or without EM contributions. But remember that it is the quantities a/y/Ni, b/y/Ni, a/\/N2 and 
bj 'y ' N2 that alter the simple quark flavor mixing angles, we should compare the mixed wavefunctions for 
77 and 77' in those fittings: 

without EM : I ) V ) 
I W) 



0.866|nn) - 0.494|ss) - 0.070|G) 
0.480|n7l) + 0.841|ss) + 0.250|G) 



(30) 



MP n Model- J to)=0.859|nfi>-0.510|a5)+0.046|G) . , , 

mr-D ivioaei. j j^) = o.503|nn) + 0.847|ss) - 0.173|G) ' (6L) 

and 

MP-C Mode,: { $ = Jj^|§ . (32, 

Generally, it shows that the glueball component in 77 is small and in 77' is relatively large. The qq 
contents are the dominant ones in their wavefunctions. In the case without EM contributions, the signs 
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for the glueball component are opposite to those with the EM contributions included. Since the glueball 
components are small, the qq mixing patterns are quite similar in these three fittings, in particular, MP-D 
and MP-C give almost the same results. We note in advance that the best fitting results are obtained in 
the MP-C model. Thus, we concentrate on the MP-C model here and try to extract information about 
the glueball candidate rj". Similar discussions can be applied to the other two schemes. 

The above three equations indicate that orthogonality between |t/) and rf) is approximately satisfied: 
(vW) = —0.6%. This would allow us to derive the mixing matrix elements for 7/" based on the unitary 
and orthogonal relation: 

/ 0.859 -0.510 0.042 \ /nn\ 
= 0.504 0.848 -0.163 ss , (33) 
V 0.090 0.140 0.986 / \G J 

where the orthogonal relation is satisfied within 5% of uncertainties. Note that an overall sign for |t/") is 
possible. 

So far, there is no much reliable information about the masses for \nn) and \ss). Also, there is no firm 
evidence for a h state (denoted as 77") as a glueball candidate and mixing with 77 and 77'. Interestingly, 
the above mixing suggests a small qq-G coupling in the pseudoscalar sector, which is consistent with 
QCD sum rule studies (45j. As a test of this mixing pattern, we can substitute the physical masses for 77, 
77' and h resonances such as 7/(1405), 7/(1475) and 7/(1835) [4(| into Eq. ([55]) to derive the pure glueball 
mass Mq, and we obtain, Mg ~ M^u due to the dominant of the glueball component in 77". Empirically, 
this allows a h glueball with much lighter masses than the prediction of LQCD [351 ] . 

The amplitude for J/ip(ip') — > tpr]" and cor)" can be expressed as 

M<fv" = 9J/il>{ii>') R [ rx 3 + (1 + r )V3 + za]F(p) , (34) 

and 

= gj/w) [(1 + 2r)x 3 + V2Rry 3 + V2z 3 ]T(p) . (35) 

For those h resonances assumed to be the 77" in Eq. l[3"3"]). predictions for their production rates in 
J/tp(ip') decays are listed in Table IVH It shows that if those states are glueball candidates, their branching 
ratios are likely at order of 10~ 3 in J ftp — > 077" and ojrj" , and at 10 -5 in ip' decays. We do not present 
the same calculations for the other two mixing schemes since they all produce similar results. 

Experimental signals for 77(1405) in J/ip radiative decays were seen at Mark III [47], |48j] and DM2 [49( ] . 
In J/tp hadronic decays DM2 reported BR < 2.5 x 10~ 4 at CL = 90% with the unknown 77(1405) ^ KKtt 
branching ratio included [H3|, and the search in Mark III I5jj gave BR{J/ip — > 077(1405) — > tftKKn) < 

1.2 x 10~ 4 {90%C.L.) and BR(J/ip -> ^7/(1405) -> ujKKtt) = (6.8 ± 1.7) x lO" 4 . 

The estimate of its total width is still controversial and ranges from tens to more than a hundred MeV 
in different decay modes 0. This at least suggests that -Bi?(?7(1405) — > KKtt) ~ 10 ~ 50%, and leads to 
BR(JN) — > 077(1405)) — » (pKKiT — < (1.49 — 7.45) x 10~ 4 . This range seems to be consistent with the 
data [Ho, [5l|. The prediction, BR(J/ip — > ^77(1405)) = 7.32 x 10~ 3 turns to be larger than 4> channel. 
A similar estimate gives, BR(J/ip — > ^77(1405) — > ujKKtt) ~ (0.73 ~ 3.66) x 10~ 3 , and can also be 
regarded as in agreement with the data from Mark III ,51]. A search for 77(1405) at BES-III in J/ip and 
tp' hadronic decays will be helpful to clarify its property. 

The same expectation can be applied to 7j(1475) and 77(1835). If they are dominated by glueball 
components, their production rate will be at order of 10 -3 as shown in Table fVTl 7/(1475), as the higher 
mass ^ in the 7/(1440) bump, tends to favor decaying into KKtt (also via ao(980)7r and K*K) [491 ]. 
As observed in experiment that the partial width for 7/(1475) — ► KKtt is about 87 MeV [2|, an estimated 
branching ratio of BR(rj(U75) -> KKtt) = 0.1 - 0.5 will lead to BR(J/ifj 07/(1475) <\>KKtt) = 
(1.33 ~ 6.65) x 10" 4 and BR{J/ip -> wt/(1475) ujKKtt) = (6.81 - 34.05) x 10~ 4 . This value is 
compatible with the production of 7/(1405), thus, should have been seen at Mark III in both u> and <fi 
channel. However, signals for 7/(1475) is only seen in <fi channel at BR(J/tp — * 07/(1475) — * (/)KKtt) < 
2.1 x 10~ 4 (90%C.L.) [fLU]. Such an observation is more consistent with the ?7(1475) being an ss-dominant 
state as the radial excited state of 77' [36l . [52| . Its vanishing branching ratio in to channel can be naturally 
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explained by the DOZI suppressions. We also note that r?(1475)'s presence in 77 — > KKtt [53j certa inly 
enhances its assignment as a radial excited ss state in analogy with 7/(1295) as the radial excited nn [36l . 

The narrow resonance X(1835) reported by BES in J/ip — + 7^(1835) — > 77r + 7r~T? / is likely to have 
J PC = h [46j], and could be the same state reported earlier in J/ip — > 7pp [HJ. Its branching ratio is 
reported to be BR(J/ip -> 7^(1835) -> 77r+7r~7/') = (2.2 ± 0.4 ± 0.4) x 10~ 4 . By assuming that it is a 
pseudoscalar glueball candidate, its production branching ratios in 077(1835) and a;7/(1835) are predicted 
in Table PvTl If 7r + n^i]' is the dominant decay channel, one would expect to have a good chance to see it 
in cl>t/(1835) channel. Since it is unlikely that a glueball state decays exclusively to ir + ir~r]' (and possibly 
pp), the predicted branching ratio BR(J/ip — > wX(1835)) = 3.66 x 10~ 3 has almost ruled out its being 
a glueball candidate though a search for the X(1835) in its recoiling u> should still be interesting. A 
number of explanations for the nature of the X(1835) were proposed in the literature. But we are not to 
go to any details here. 

With the above experimental observation, our results turn to favor the 7/(1405) being a pseudoscalar 
glueball candidate. 



3. Mixing in an old perturbation theory 

In this scheme, the masses of the nn, ss and glueball G are treated as parameters along with the 
"perturbative" transition amplitudes f a and fb- What we refer to as "perturbative" here is that both 
f a and fb have values much smaller than the mass differences between those mixed states. The fitted 
results in Table fVl indeed satisfy this requirement. Parameters f a is fitted to be around 76 MeV, which 
suggests a rather small glueball component in 77 and 77' wavefunctions. In contrast, the Tin and ss mixing 
strength is slightly larger, i.e. f b ~ 94 MeV. The fitted masses M nn ~ 0.658 GeV and M sS ~ 0.853 GeV 
are located between the 77 and 7/', which is consistent with the expectation of the quark- flavor mixing 
picture with a mixing ang le (j> = fl + arctanv^ ~ 34.7°. The fitted mass for the glueball is about 1.4 GeV, 
which is determined by the assumption that the 7/(1405) is a pseudoscalar glueball candidate. Since the 
transition amplitude f a is relatively small, the mixing does not bring significant differences between the 
pure and physical glueball masses. 

The wavefunctions for 77, 77' and 77" in the three different considerations of the EM transitions are 

(36) 



(37) 





\v) - 


= 0.864|n7i) - 


- 0.479|ss) 


- 0.157|G) 


without EM : j 


W) - 


= 0.480|nn) - 


f 0.864 ss) 


-0.152 G) 




, IV) = 


= 0.176 nn) - 


f 0.170 ss) 


+ 0.970 G) 




f 1") : 


= 0.864|nn) - 


- 0.478| ss) 


-0.155|G) 


MP-D Model: j 


W) 


= 0.479 nn) - 


f 0.865 ss) 


- 0.150 G) 




[ W) 


= 0.174 nn) - 


f 0.168 ss) 


+ 0.970 G) 




{ \v) 


= 0.868|nn) 


- 0.473|ss) 


-0.153|G) 


MP-C Model: < 


W) 


= 0.473 nn) 


+- 0.869 ss) 


- 0.147 G) 




{ W) 


= 0.171 7in) 


+ 0.165 ss) 


+ 0.971 G) 



and 

f 1 77) = 0.868|nn) - 0.473 |ss) - 0.153|G) 

(38) 

In comparison with the first two mixing schemes, this approach in the framework of old perturbation 
theory produces a similar mixing pattern as that in Scheme-II except that the glueball component in the 
77 wavefunction is quite significant, e.g. as shown in Eq. 



4- The branching ratios for J/ip(ip') —> VP and violation of the "12% rule" 

The fitted branching ratios for J/ip(ip') — > VP are listed in Tables IVTT1 and IVIII1 and the results for the 
isospin violated channels arc also included as a comparison. It shows that all these three parametrizations 
can reproduce the data quite well though there are different features arising from the fitted results. 
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One predominant feature is that in all the schemes, parameter gj/^r^n is found to have overall consistent 
values for both J/tp and tp'. As pointed out earlier that the relatively small value of gj/^w) for tp' suggests 
an overall suppression of the tp' — > 3g form factor. With the destructive interferences from the relatively 
large EM contributions, the branching ratios for tp' — > pir are further suppressed and this leads to the 
abnormally small branching ratio fractions between the exclusive decays of J/tp and tp' . Numerically, 
this explains why the 12% rule is violated in J/tp(tp') — > pir. Nonetheless, such a mechanism is rather 
independent of the final state hadron, thus, should be more generally recognized in other exclusive 
channels. This turns to be true. For instance, the large branching ratio difference between the charged 
and neutral K*K + c.c. channels highlights the interferences from the EM transitions where the 
relative phases to the strong amplitudes are consistent with the expectations for the pn channel (2pj . 

In line with the overall good agreement of the fitting results to the data is an apparent deviation in 
tp' — ► lot]' in Scheme-II and III. The numerical fitting in Scheme-II gives BR(tp' — > ujrf) — 9.49 x 10~ 8 

+ 2 5 r-i 

in contrast with the data (3.2 ) x 10 -5 [2|. The significant deviation from the experimental central 

value is allowed by the associated large errors. In fact, this channel bears almost the largest uncertainties 
in the datum set. The small values from the numerical fitting also reflect the importance of the EM 
interferences. Note that in the fitting with only strong transitions, the branching ratio for lot]' has 
already turned to be smaller than the data. With the EM contributions, which are likely to interfere 
destructively, this channel is further suppressed. 

In contrast, the CKM mixing scheme is able to reproduce the tp' — > lot]' branching ratio. This is because 
the qq and G components are in phase in the rj' wavefunctions in Eqs. (|2"T|l - P§|) . Since the tp' — > lot]' 
decay shows large sensitivities to the mixing scheme, it is extremely interesting to have more precise data 
for this channel as a directly test of the mixing schemes proposed here. As BES-II may not be able to 
do any better on this than been published [231] . CLEO-c with a newly taken 25 million tp' events can 
presumably clarify this [55| . 

In Table llXj we present the branching ratio fractions R for all the exclusive decay channels for the MP-C 
model. The extracted ratios from experimental data are also listed. Apart from those three channels, cor)', 
prf and (pir, of which the data still have large uncertainties, the overall agreement is actually quite well. It 
clearly shows that the "12% rule" is badly violated in those exclusive decay channels, and the transition 
amplitudes are no longer under control of pQCD leading twist [l|. The power suppression due to the 
violation of the hadronic helicity conservation in pQCD will also be contaminated by other processes 
which are much non-perturbative, hence the pQCD-expected simple rule does not hold anymore. Note 
that it is only for those isospin violated channels with exclusive EM transition as leading contributions, 
may this simple rule be partly retained [l3[ . 



IV. SUMMARY 



In this work, we revisit J/tp(tp') — > VP in a parametrization model for the charmonium strong decays 
and a VMD model for the EM decay contributions. By explicitly defining the SU(3) flavor symmetry 
breaking and DOZI violation parameters, we obtain an overall good description of the present available 
data. It shows that a reliable calculation of the EM contributions is important for understanding the 
overall suppression of the tp' — > 3g form factors. Our calculations suggest that pn channel is not very much 
abnormal compared to other VP channels, and similar phenomena appear in K*K as well. Meanwhile, 
we strongly urge an improved experimental measurement of the tp' — > tor]' as an additional evidence for 
the EM interferences. Although it is not for this analysis to answer why tp' — > 3g is strongly suppressed, 
our results identify the roles played by the strong and EM transitions in J/tp(tp') — > VP, and provide 
some insights into the long-standing "pir puzzle" . 

Since the EM contributions are independently constrained by the available experimental data [2j, the 
parameters determined for the pseudoscalar in the numerical study, in turn, can be examined by those 
data. In particular, we find that r\ and rf allow a small glue component, which can be referred to the 
higher Fock state contributions due to the U(1)a anomaly. This gives rise to the correlated scenario 
between the octet-singlet mixing angle and the decay constants as addressed by Feldmann et al [HI, H3] ■ 

We are also interested in the possibility of a higher glue-dominant state as a h glueball candidate. 
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Indeed, based on the fact that only a comparatively small glueball component exists in rj and rf, we 
find that a ^ glueball which mixes with 77 and rj' is likely to have nearly pure glueball configuration. 
Although the obtained mixing matrix cannot pin down the mass for a glueball state, the glueball dom- 
inance suggests that such a glueball candidate will have large production rate in both <frq" and corj" at 
10 -3 . This enhance the assignment that ^(1405) is the h glueball candidate if no signals for ry(1475) 
and 77(1835) appear simultaneously in their productions with <fr and lj in J/ip(ip') decays. High-statistics 
search for their signals in <f> and u> channel at BESIII will be able to clarify these results. 
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Coupling const. 


e/fv 


Values (xlO -2 ) 


Total width of V 


BR(V -> e+O 


e/fp 




6.05 


146.4 MeV 


(4.70 ± 0.08) x 10" 5 


e/f u 




1.78 


8.49 MeV 


(7.18 ±0.12) x 10~ 5 


e/U 




2.26 


4.26 MeV 


(2.97 ±0.04) x 10~ 4 


e /fj/i> 




2.71 


93.4 keV 


(5.94 ± 0.06)% 


e/U> 




1.65 


337 keV 


(7.35 ±0.18) x 10~ 3 



TABLE I: The coupling constant e/fv determined in V 
PDG2006 0]. 



The data for branching ratios are from 



Decay channels 


Transition amplitude M = {Mi +M2+M3) 


H 


gj/ il (i,>)R[V2rx 1 ± R(l + r)yi + Zi]F(P) 


<H 


gj /4 ,(4,>)R[V2rx2 + R(l + r)y 2 + Za]^(P) 




gj/iPO-nK 1 + 2r 0^i + V2Rr yi + </2zi]F(P) 


UJTj' 


ffj/WjK 1 + 2r ) X2 + V^Rryi + V2z 2 ]F(P) 


A° 




p + ir~ or p~ir + 


gj/i,a,')^(P) 


K*°K° or K*°K° 


g.j/i,(4-')RH p ) 


K*+K~ or K*~K+ 





TABLE II: General expressions for the transition amplitudes for J/ip(tp') — * VP via strong interactions. 



15 





without EM 


MP-D 


MP-C 






J/i> 




J/iP 


i>' 


r 


-0.265 ±0.023 


-0.185 ±0.194 


-0.270 ± 0.024 


-0.080 ±0.155 


-0.273 ±0.024 


-0.054 ± 0.162 


R 


0.660 ± 0.023 


1.400 ±0.112 


0.667 ± 0.024 


1.293 ±0.163 


0.667 ±0.024 


1.233 ±0.238 


5j/w)( x10 ~ 3 ) 


18.94 ±0.64 


1.66 ±0.15 


18.62 ±0.63 


2.06 ±0.25 


18.66 ±0.63 


2.22 ± 0.49 


s 






75.0° ± 5.2° 


143.9° ± 14.6° 


73.7° ± 5.8° 


159.5° ±29.9° 


012 


156.2° ± 2.0° 


154.7° ± 2.4° 


154.3° ± 2.4° 


013 


179.8° ± 0.7° 


180.0° ± 0.3° 


180.0° ±0.2° 


023 


9.4° ± 1.7° 


9.5° ± 1.9° 


9.7° ± 1.9° 


X 2 /d.o.f 


37.0/9 


9.8/11 


9.1/11 



TABLE III: Parameters introduced in the overall fitting of J/ip(tp') — > VP in Scheme-I, i.e. the CKM approach. 





without EM 


MP-D 


MP-C 


J/V 


</>' 


J/V 




J/V 


v/ 


r 


-0.183 ±0.034 


0.101 ±0.237 


-0.308 ±0.037 


-0.237 ±0.155 


-0.307 ±0.038 


-0.206 ±0.129 


R 


0.661 ±0.024 


1.400 ±0.106 


0.672 ± 0.027 


1.098 ±0.232 


0.674 ± 0.028 


1.065 ±0.181 


5j/VW')( x10_3 ) 


18.91 ±0.64 


1.63 ±0.15 


18.45 ±0.69 


2.40 ±0.53 


18.45 ±0.70 


2.54 ±0.44 


8 






74.0° ± 5.2° 


30.7° ± 20.4° 


72.3° ± 6.0° 


11.0° ±61.8° 





-25.0° ±3.7° 


-24.0° ±0.8° 


-24.0° ±0.9° 


a 


0.673 ±0.159 


(3.15 ±0.90) x 10~ 2 


(3.44 ± 1.03) x 10~ 2 


b 


-0.127 ±0.013 


(4.03 ±0.60) x 10~ 3 


(4.13 ±0.65) x 10~ 3 


xVd-o.f 


41.1/9 


11.2/11 


10.7/11 



TABLE IV: Parameters introduced in the overall fitting of J/ip(ip') — > VP in Scheme-II, i.e. qq-G mixing due to 
higher Fock state contributions. 





without EM 


MP-D 


MP-C 


J/.A 


v/ 


JH> 


i>' 


J/V 


</>' 


r 


-0.035 ± 0.056 


0.115 ±0.193 


-0.037 ±0.055 


0.031 ±0.151 


-0.044 ±0.057 


0.041 ±0.141 


R 


0.698 ± 0.025 


1.400 ±0.297 


0.715 ±0.027 


1.079 ±0.155 


0.718 ±0.028 


1.047 ±0.158 


3j/<W)( xlCr3 ) 


18.05 ±0.64 


1.65 ±0.15 


17.53 ± 0.64 


2.44 ± 0.34 


17.52 ±0.66 


2.57 ±0.42 


8 






74.9° ± 5.07° 


149.2° ± 19.6° 


73.8° ± 5.8° 


168.2° ± 136.7° 


M nn (GeV) 


0.658 ±0.083 


0.658 ± 0.055 


0.657 ±0.081 


Ms- S (GeV) 


0.855 ± 0.076 


0.853 ± 0.057 


0.853 ± 0.077 


M G (GeV) 


1.389 ±0.093 


1.390 ±0.093 


1.391 ±0.095 


fa (GeV) 


0.077 ± 0.039 


0.076 ± 0.022 


0.076 ± 0.038 


U (GeV) 


0.094 ±0.033 


0.093 ± 0.034 


0.091 ±0.036 


X 2 /d.o.f 


46.3/13 


18.9/13 


18.3/13 



TABLE V: Parameters introduced in the overall fitting of J/tp(xp') — » VP in Scheme-Ill, i.e. mixing via an old 
perturbation theory. 
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v" 


BR(J/^ -> Vr)") (xKT 3 ) 


BR{4,' -» V-q") (xlQ- 5 ) 


H" 


uir] 


H' 


// 

OJT] 


77(1405) 


1.49 


7.32 


4.07 


7.17 


77(1475) 


1.33 


6.81 


3.96 


7.03 


77(1835) 


0.44 


3.66 


3.09 


5.93 



TABLE VI: Predictions for production of 77" in J/ip(tp') -> VP by assuming it is 77(1405), 77(1475) and 77(1835), 
respectively, in Scheme-II. Calculations by Scheme-I and III give similar results with the same order of magnitude. 



Decay channels 


without EM 


Scheme-I 
(MP-C) 


Scheme-II 
(MP-D) 


Scheme-II 
(MP-C) 


Scheme-Ill 
(MP-C) 


Exp. data 


„o 

p 7T 


5.64 


X 


10" 


3 


5.76 x 10" 3 


5.64 x 10 -3 


5.65 x 


10" 3 


5.09 x 


10" 


3 


(5.6 ±0.7) x 10~ 3 


pTT 


1.69 


X 


10" 


2 


1.73 x 10" 2 


1.70 x 10~ 2 


1.69 x 


10 -2 


1.53 x 


lo- 


2 


(1.69 ±0.15) x 10" 2 


tor] 


1.55 


X 


10" 


3 


1.53 x 10- 3 


1.60 x 10 -3 


1.60 x 


10 -3 


1.75 x 


ur 


3 


(1.74 ±0.20) x 10" 3 


LOT]' 


1.82 


X 


10" 


4 


1.82 x 10~ 4 


1.83 x 10~ 4 


1.83 x 


10~ 4 


1.80 x 


10" 


4 


(1.82 ±0.21) x 10~ 4 




6.45 


X 


10" 


4 


6.40 x 10~ 4 


6.91 x 10~ 4 


6.94 x 


10" 4 


5.92 x 


10" 


4 


(7.4 ±0.8) x 10" 4 




4.00 


X 


10" 


4 


4.22 x 10" 4 


4.00 x 10~ 4 


4.00 x 


10" 4 


3.17 x 


10" 


4 


(4.0 ±0.7) x 10" 4 


K*+K~ +c.c. 


4.67 


X 


10" 


3 


5.01 x 10~ 3 


5.03 x 10~ 3 


5.03 x 


10 -3 


5.11 x 


10" 


3 


(5.0 ±0.4) x 10" 3 


K*°K° + c.c. 


4.66 


X 


10" 


3 


4.30 x 10 -3 


4.25 x 10 -3 


4.23 x 


10 -3 


4.39 x 


10" 


3 


(4.2 ±0.4) x 10 -3 


PV 










1.5 x 10 -4 


1.1 x 10 -4 




1.5 x 10~ 4 






(1.93 ±0.23) x 10 -4 


pv' 










7.9 x 10~ 5 


4.2 x 10 -5 




7.9 x 10" 5 






(1.05 ±0.18) x 10~ 4 












3.3 x 10 -4 


4.2 x 10 -4 




3.3 x 10 -4 






(4.5 ±0.5) x 10~ 4 


<f>ir 










9.9 x 10 -7 


8.0 x 10" 7 




9.9 x 10 -7 






< 6.4 x 10 -6 



TABLE VII: Fitted branching ratios for J/tp — > VP in different parametrization schemes for the isoscalar mixings. 
Column (without EM) is for results without EM contributions; Columns of MP-C correspond to processes Fig.[2ja) 
and (b) in a constructive phase with an effective mass A = 0.616 GeV; and Column MP-D corresponds to (a) and 
(b) in a destructive phase with A = 0.65 GeV, which are the same as in Ref. [ljj. The isospin violated channels 
(last four channels) are also listed. The experimental branching ratios are from PDG2006 Q{. 
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Decay channels 


without EM 


Scheme-I 


Scheme-II 


Scheme-II 


Scheme-Ill 


Exp. data 












(MP-C) 




(MP-D) 


(^lvlr^-U ) 


(MP-C) 







p IT 


0.89 


X 


10" 


5 


0.71 x 10" 


-5 


0.88 x 


10~ 5 


0.97 x 10" 5 


1.02 x 10" 


5 




pir 


2.68 


X 


10" 


5 


2.03 x 10 


-5 


2.58 x 


10~ 5 


2.78 x 10 -5 


2.92 x 10" 


5 


(3.2 ± 1.2) x 10~ 5 


Lurj 


6.36 


X 


10" 


6 


6.08 x 10" 


-6 


5.63 x 


10" 6 


5.35 x 10 -6 


4.97 x 10" 


6 


< 1.1 x 10 -5 


LJT]' 


1.35 


X 


10" 


6 


9.68 x 10 


-6 


1.96 x 


io- 7 


9.49 x 10 -8 


5.21 x 10" 


7 


_1_9 K 

(3.2 +2 " 5 ) x nn 5 

-2.1 


<Pv 


u.yo 


X 




5 


1 79 Y 1(1" 


-5 


2.02 x 


10~ 5 


2.04 x 10" 5 


1 Of! y 1(1" 


5 


_i_i n 

(2.8 + ) x IO" 5 
-0.8 


H 


1.75 


X 


10" 


5 


2.11 x 10 


-5 


2.13 x 


lO^ 5 


2.41 x 10~ 5 


2.80 x 10" 


5 


(3.1 ±1.6) x 10~ 5 


K*+K~ +c.c. 


3.50 


X 


10" 


5 


2.28 x 10 


-5 


2.07 x 


10" 5 


2.01 x 10~ 5 


2.00 x 10" 


5 


(1.7 +0 " 8 ) x IO" 5 
—0.7 


K* a K° + c.c. 


3.50 


X 


10" 


5 


1.18 x 10 


-4 


1.18 x 


IO" 4 


1.17 x 10~ 4 


1.13 x 10" 


4 


(1.09 ±0.20) x 10~ 4 


PV 










1.4 x 10" 


5 


9.4 x 10 -6 


1.4 x 10 -5 




(2.2 ± 0.6) x IO -6 


prj' 










7.4 x 10" 


6 


3.9 x 10 -6 


7.4 x 10 -6 




(1.9 +L7 ) x IO" 5 
-1.2 


UJTV 










3.0 x 10" 


5 


3.9 x 10 -5 


3.0 x 10~ 5 




(2.1 ±0.6) x IO -5 


(f>1T 










7.3 x 10" 


8 


9.6 x 10~ 8 


7.3 x 10" 8 




< 4 x 10~ 6 



TABLE VIII: Fitted branching ratios for ip' -v VP. The notations are the same as Table IVIII The stars 
in p°7r° channel denotes the unavailability of the data. 



Decay channels 


Scheme-I (%) 


Scheme-II (%) 


Scheme-Ill (%) 


Exp. data (%) 


pn 


0.12 


0.15 


0.19 


0.2 ±0.1 


LOT] 


0.40 


0.35 


0.28 


< 0.6 ±0.1 


LOTj 


5.33 


0.11 


0.29 


18.5 ± 13.2 


<f>V 


2.78 


2.93 


3.30 


4.1 ± 1.6 




5.00 


5.34 


8.86 


8.7 ±5.5 


K *+ K - +cx . 


0.45 


0.41 


0.39 


0.4 ± 0.2 


K*°K° + c.c. 


2.74 


2.79 


2.67 


2.7 ±0.7 


PV 




8.97 




11.5 ±5.0 


prj' 




9.44 




23.5 ± 17.8 


um 




9.01 




5.0 ± 1.8 


<f)1T 




7.41 




< 62.5 



TABLE IX: Branching ratio fractions for all VP channels in the MP-C model, 
also included. 



The isospin violated channels are 
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FIG. 1: Schematic diagrams for J/ip — > <j>P via strong interaction, where the production of different components 
of the pseudoscalar P is demonstrated via (a): SOZI process; (b) DOZI process; and (c) glueball production. 
Similar prcesses apply to other VP channels as described in the text. 
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FIG. 2: Schematic diagrams for J/ip(tp') — > 7* — > VP. 



